function b = newton(a, x, y)
	 % Intrari:
	 % a = abscisa in care se calculeaza polinomul
	 % x = tabloul absciselor celor n puncte
	 % y = tabloul ordonatelor celor n puncte
         % Iesire:
         % b-valoarea polinomului de interpolare in a

	 n=length(x);
	 c=difdiv(x,y);
	 b=c(1);
	 p=1;
	for i=2:n
		p=p*(a-x(i-1));
		b=b+p*c(i);
	endfor

endfunction

